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NOMENCLATURE 



Symbo 1 
b 

b l 

b 2 

B A (b 1 + b 2 ) 

B D ( b l +b 2 ^ 
B E (b) 

B E ( b i ,b 2 ) 

b| ( b 1 ,b 2 ) 

E 

E 0 

E i 

E n ( x ) 

F 

I (b 1 ,E ,fi) 
I S (b 1 ,E ,3) 



Descri pti on 

Thickness in mean free paths 
(b = px ) 

Thickness of the first layer in mean 
free paths 

Thickness of the second layer in mean 
free paths 

The two-medium e ne rgy- ab so rp t i on 
buildup factor 

The two-medium dose buildup factor 

The single- medium energy buildup 
factor 

The energy buildup factor at the 
interface between the two layers 

The scattered energy buildup factor 

Photon energy 

The initial photon energy 

The average photon energy in the 
i th group 

The exponential integral of order n 

The fractional deviation 

The angular energy flux density at the 
interface between the slabs, defined 
as E . (j, ( b-j ,E,3) c 



Dimension 

None 

None 

None 

None 

None 
N one 

None 

None 

MeV 

MeV 

MeV 

None 

None 

MeV 

^-sec-MeV-STER 



The scattered angular energy flux Me V 

at the interface between the slabs cm 2 _ sec _^ e \/-STER 
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Synbo 1 Descri ption 

^(b^b^) The scattered energy flux at the 
interface between the slabs 



Dimension 

MeV 

2 

cm -sec 



I s (b 1 +b 2 ) 



The scattered energy flux at the 
outer edge of the second layer 



MeV 



2 

cm 



-sec 



I (bi+bp^E) The scattered energy flux spectrum 
at the outer edge of the second 
1 aye r 

I (bpb^E) The ^egendre^ expansion coefficients 



I u (b n + b 2 ) 



of I ( b i , E , w) 

The unscattered energy flux at the 
outer edge of the second layer 



I SS (b^+b 2 ) The energy flux which scatters in 
both 1 aye r s 



s u * 

I ( b" i + b 2 ) The energy flux which scatters in 
the first layer but not in the 
second 1 ayer 



u s 

I ( b -| + b ^ ) The energy flux which does not 
scatter in the first layer but 
scatters in the second layer 



I?( b l ) 



ij:(b 1 ,b 2 ) 



MeV 



cm -sec-MeV 
MeV 

2 

cm -sec-MeV 

MeV 

2 

cm -sec 

MeV 

2 

cm -sec 
MeV 



The scattered energy flux which 
results at the interface when the 
second layer is replaced by a vacuum 

The scattered energy flux at the 
interface which is reflected from 
the second 1 ayer 



MeV 



cm -sec 



MeV 



cm -sec 



MeV 



cm -sec 



s u 

I n ( b i + b 2 ) The energy flux which scatters in the None 
first layer but not the second layer, 
normalized to an interface flux of 
lMeV/cm-sec. 



M (b) 

£ 



ri^ U (b 1 +b 2 ) = I SU (b 1 +b 2 )/I S (b 1 ,b 2 )] 

Defined by equation (3.12) None 

The rest mass of an electron gm 
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Symbol 


Description 


Dimension 


N 


The atom density 


a toms /cm'* 


V"> 


Legendre polynomial of order £ 


None 


R(b 1 ,b 2 ) 


The ratio of scattered to unsca£- u 
tered energy flux albedo (R = a^/a^) 


None 


- V 

r 


An independent variable for position 


cm 


r 

e 


Ttie classical radius of an electron 


cm 


s 


The standard deviation 


None 


S N 


The order of angular quadrature 


N one 


w 


Cosine of the angle of emergence 
(see Figure 3.1) 


Radians 


X 


Distance in slab geometry 


cm 


X 1 


Thickness of the first layer 


cm 


X 2 


Thickness of the second layer 


cm 


z 


Atomic number 


None 


Z 1 


Atomic number of the material in the 
first 1 aye r 


None 


z 2 


Atomic number of the material in the 
second layer 


None 


a^(b) 


The unscattered energy flux albedo 


None 


aj( b i .b 2 ) 


The scattered energy flux albedo 


None 


6 


The Dirac delta function 


None 


ae. 


Width of the ith energy group 


MeV 


X 


Photon wavelength in Compton units 


. 02426° A 


P 


The linear attenuation coefficient 
(p = No) 


-1 

cm 


y l 


The linear attenuation coefficient 
for the material in the first layer 


cm * 


u 2 


The linear attenuation coefficient for 
the material in the second layer 


cm * 



x x i i i 



Symbol 

u a 

^en 



Description 

The linear absorption coefficient 

The linear energy-absorption 
coefficient 



D i m e n s i o n 



cm 



-1 



cm 



-1 



v en /0 



The mass energy-absorption coef- 
ficient 



c m / g m 



The microscopic photoelectric and 
pair production absorption cross 
section 



cm 



°t 

h-g 

s a 



The microscopic total cross section 
(total = absorption + scattering) 

The P ^ down scattering cross sections 



cm 



cm 



-1 



E ( r , E ' *-E , ft 1 ) The differential cross sectio^n cm 

for scattering from the direction u' MeV-STER 
to ft and from the energy E' to E per 
unit solid angle and energy range 

I (r,E'->-E) The Legendre expansion coefficient of cm 

s Z s (f,E'-E,sS'^) Mi\T 

<|> ( b, , E , /) The angular photon flux, defined so Photons 
that 4>(b-, ,E,T<:)dE dw gives the number 
of photons at b ^ , with energy E in the 
range dE, goingin the direction ft 
within the element of solid angle d , 
which cross in unit time a unit dif- 
ferential element of area whose normal 
is in the direction 



cm 



- se c -Me V ■ 
STER. 
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ABSTRACT 

Mam'scalco, James Andrew. Ph.D., Purdue University, 

May 1973. AN ANALYTICAL TREATMENT OF PHOTON TRANSPORT IN 
TWO-LAYERED MEDIA. Major Professor: H. E. Hungerford. 

In this thesis the transport of photons in two-layered 
slabs is examined, and then evaluated by analytically ex- 
pressing the two-medium buildup factor in terms of single 
medium buildup factors and parameters which in turn can be 
characterized as functions of atomic number. The word 
analytical is used to indicate that the parameters in this 
formulation are not empirical but rather they are directly 
related to the basic physics of transport phenomena in two - 
1 ayered slabs. 

In applying the two-medium buildup factor formulation, 
we assume that the single medium buildup factors and albedos 
are available and then seek to determine the two-medium 
parameters. These quantities have been evaluated by using a 
one-dimens i ona 1 discrete ordinates transport code (ANISN) to 
calculate the attenuation of photons in two-layered slabs 
of water, iron, and lead. Independent verification of some 
of the ANISN results was made using the stochastic transport 
code PUGTII. The materials water, iron, and lead were chosen 
because they are commonly used together in pair combinations 
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and because they represent the light, intermediate, and 
heavy elements respectively. 

In the process of tabulating and plotting the results, 
it was observed that the two-medium parameters are slowly 
varying functions of the first medium and smoothly varying 
functions of the second medium; moreover, a detailed examin- 
ation of the plotted results (Appendix D) reveals that the 
curves for each of the two-medium parameters have the same 
shape (i.e., functional dependence) for different first- 
layer thicknesses and/or materials. Because of this func- 
tional behavior, these two-medium parameters have been 
accurately extrapolated to a large number of two-material 
combinations by simple linear interpolation on atomic num- 
ber . 

Energy and dose buildup factors for more than 15 differ- 
ent two-material combinations have been calculated with the 
two-medium buildup factor formulation and compared to cor- 
responding ANISN results. These comparisons clearly indi- 
cate that the errors resulting from obtaining buildup factors 
with this formulation are less than 5%. 
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1. INTRODUCTION 

1.1. Statement of the Problem 
The advent of the nuclear age has resulted in the de- 
velopment of sophisticated nuclear devices and systems in 
two vital areas of our material well-being: (1) national 

security; and (2) energy. The weapons and nuclear power 
reactors conceived and developed in the many military and 
civilian programs over the past quarter of a century have 
been possible only because of our advancement in the under- 
standing of nuclear radiation transport processes within ma- 
terial media. As part of this theoretical undertaking, a 
considerable amount of effort has been devoted to developing 
methods of calculating the deep penetrations of photons in 
material media. These radiation transport calculations 
usually involve the solution of the Boltzmann transport equa- 
tion in one of its various forms. 

The Boltzmann transport equation is a statistical- 
mechanical statement of the balance and flow of radiation in 
terms of position, energy, direction, and time. As such, it 
is an integro-differential equation which can not be solved 
exactly for any but a few highly restricted cases; however, 
numerous methods have been developed and used to provide 
approximate solutions. They include techniques utilizing: 
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1. Successive integrations over trial solutions, 

2. Polynomial expansions such as spherical harmonics, 

3. Fourier transforms, 

4. Diffusion theory, 

5. Variational techniques and perturbation theory, 

6. Expansions in terms of moments, 

7 . Invariant imbedding, 

8. Discrete ordinates S fJ theory, and 

9. Stochastic Solutions (Lattice Transport and Monte 
Carlo). 

Some of these methods most applicable to Dhoton transport 
(the topic of interest in this investigation) will be described 
in the next chapter. 

It has become customary to present the results of photon 
transport theory calculations in the form of a buildup fac- 
tor. The buildup factor is defined as the ratio of the de- 
sired quantity which is a c ha r ac te r i s t i c of the total gamma- 
ray flux density to a similar quantity which is a character- 
istic of the unscattered flux density. For example, the 
energy buildup factor is defined as the ratio of the total 
energy flux (direct plus scattered) to the unscattered energy 
flux (direct). The principal advantage of expressing the 
results of photon attenuation calculations in the form of 
buildup factors is that these quantities are rather slowly 
varying functions of attenuation distance, photon source 
energy, and average atomic composition. 
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The penetration of gamma-rays in an infinite, homogeneous 

medium has been explored in detail using the moments meth- 
T 1 2341 

0 d , ’ ’ 5 J a technique particularly suited to photon trans- 

port in infinite media. Goldstein and Wilkins^ 1 ! have used 
this method to compile an extensive collection of buildup 
factor data which for 20 years has been and is still con- 
sidered to be the primary source of buildup factor data for 
point isotropic and plane monodirectional sources in infinite 
media. Of more practical interest to shielders is the trans- 
mission of photons through finite slabs of material. Much 
work has been done on single (one-layer) slabs, C5,6,7,8] - n _ 
eluding calculations of reflection coefficients and transmis- 
sion buildup factors for several different source distribu- 
tions (isotropic, perpendicular, and oblique incidence). 

A small amount of work has been done on two-layer 
slabs? 9 ’ 10 ’ 11 ’ 12 ! However the object of such studies has 
been to derive corrections to single medium buildup factors 
and/or empirical formulas^ 12 ’ 1 ^’ 1 ^’ 1 ^ ^ which work only in a 
given set of circumstances. No investigation that the author 
is aware of has focused on expressing the two-medium buildup 
factor in terms of general parameters which are directly re- 
lated to the basic physics of transport phenomena in such 

1 aye r s . 

The objective of this thesis is to investigate the trans- 
port of gamma-rays in two-layer slabs and to develop an 
analytical expression for two-medium buildup factors in terms 
of single-medium buildup factors and parameters which, in 
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turn, can be characterized as a function of atomic number. 

The reason for expressing the two-layer results in terms of 
one-layer results is that one-layer results are readily 
available in the open literature, and two-layer buildup 
factors so expressed will be of great practical benefit. 

The present theoretical investigation of photon transport 
in two-layered slabs was conducted by first thoroughly re- 
viewing the literature for all existent one- and two-layer 
buildup and transmission information. As a result of this 
extensive search, it became evident that the problem of two- 
layer buildup factors could best be handled by utilizing 
both discrete ordinates and stochastic transport techniques 
to calculate, as accurately as possible, the attenuation of 
gamma-rays in two-layered material combinations of water, 
iron, and lead. These elements were chosen for the study 
because they are commonly used together in pair combinations, 
and they represent the light, intermediate, and heavy elements 
respectively. The results of these calculations were then 
used to determine a set of parameters required in the formu- 
lation of the two-material buildup factors. 

As has already been pointed out, a number of computa- 
tional techniques are available for determining the transport 
of gamma-rays in two-layered media. However, these tech- 
niques usually involve the extensive use of a large computing 
facility, and when a large number of calculations on complex 
shielding configurations are to be undertaken, this method 
is unduly slow, awkward, and costly. We looked, therefore, 
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for ways to simplify the approach for others, performing 
in advance (as it were) the long computer calculations for 
them. The unique feature of the analytical formulation 
derived in this work is that two-material buildup factors 
can be easily hand-calculated, using fairly simple formulas 
for most of the two-material combinations of interest to 
shield designers. In fact, the major limitation on the 
number of two-material combinations which can be investigated 
is the lack of one- medium data. 
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2. THE TRANSPORT OF GAMMA RAYS: 

REVIEW OF THE LITERATURE AND METHODS 

The attenuation of gamma-rays in material media is rela- 
tively simple to predict when each interaction process 
(collision) results in the disappearance of a photon, i.e., 
the medium is purely absorptive. For example, consider the 
case of a monodirectional beam of gamma -rays of energy E in 
such an absorptive medium; then, the change in the flux in- 
tensity of the beam dl(x), as it traverses the medium, is 
given by 



dx 



P I (x ) 



cm 



MeV 

3 



-sec 



( 2.1 ) 



where the constant of proportionality p is known as the 
linear attenuation coefficient. The solution of equation 
(2.1) is the well-known Lambert's law of absorption, 



I (x ) = I e _VJX - i r - V - ■ . (2.2) 

cm - sec 

Unfortunately, Lambert's Law is not an accurate description 
of photon attenuation in material media. First, the law 
assumes that the medium is purely absorptive and does not 
account for scattering. Second, the law does not account for 
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the energy change which results from photons of lower energy 
contributing to the flux at x. 

Lambert's law can be corrected to account for scattered 
radiation by using a multiplicative factor called the 
buildup factor. The term "buildup factor" was first intro- 
duced during the Manhattan Project. Its usefulness stems 
from the fact that radiation quantities resulting from un- 
scattered photons, i.e., photons which emerge without suf- 
fering any collisions, are relatively simple to calculate. 
There are many different quantities of interest in shield 
design (particle flux, energy flux, dose rate, etc.) and 
we can define a buildup factor for each of them. To show 
how this can be done, let the superscripts u and s refer to 
unscattered and scattered photons respectively and let the 
quantity of interest be represented by the functional F(I); 
then, the buildup factor for the quantity F(I) is defined 
bytn 



F ( I } = B f F(I U ) . (2.3) 

In words, equation (2.3) defines the buildup factor as the 
ratio of a desired quantity which is a characteri sti c of the 
total flux density (unscattered plus scattered) to a similar 
quantity which is a characteristic of the unscattered flux 
density. From this definition, it is obvious that the build- 
up factor concept is meaningless when the unscattered quan- 
tity does not exist. 
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Buildup factors of interest in shield design include: 

1. The energy buildup factor B^. (with F(I) represent- 
ing the integrated energy flux) defined as 



B E ( x ) : 


Jr I(x,E) dE 

: H • (2.4) 

J E I ( x , E ) dE 


The total energy 


flux I ( x , E ) is the sum of the scattered 



and unscattered contributions. That is, 



I ( x , E ) 


= I U ( x , E ) + I S ( x , E ) (2.5) 



Hence , 



B E ( x ) = 


Jr I ( x , E ) dE 

; 1 + ~ r • (2.6) 

/ E I ( x , E ) dE 



2. The dose buildup factor Bp (with F(I) reoresenting 
the dose-rate in air) defined as 



B d (x) = 


/ U y (E) I s (x,E) dE 

: 1 + a i r U ’ ( 2 • 7 ) 

J E V en (E) I U ( x , E ) dE 



Q 1 p 

where P en (E) is the energy absorption coefficient of air. 

3. The energy absorption buildup factor B ^ (with F(I) 



representing the 


energy density absorbed by the medium 


from the photons 


in unit time) defined as 


CD 

X 


/ r u,„ (E) I S (x,E) dE 
- 1 + E e " „ , (2.8) 
/ E u en (E) I U ( x , E ) dE 
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where (E) is the energy absorption coefficient for the 

medium. It should be mentioned that in all of the defini- 
tions cited, the functional F depends only on the energy 
flux; therefore, for a monoenergeti c source, F(I U ) is given 
by a constant, F , times the integrated, unscattered energy 
flux I u (x) . In these cases, the buildup factor is simply 
defined by 

F ( I ) - B p F u I u (x) • (2.9) 



2.1. Transport Theory 

The fundamental equation governing the transport of 
gamma rays in material media is an i n t egro - d i f f eren t i a 1 
equation known as the Boltzmann transport equation. It is 
basically a bookkeeping statement which accounts for photon 
increases and decreases in a given increment of space, 
angle, energy, and time. For shield analyses, a steady 
state is normally assumed and the time variable is not re- 
quired. Briefly, the equation of continuity is applied to a 
s i x- d i mens i ona 1 phase space consisting of three spatial 
coordinates r, two directional coordinates Q and one energy 
coordinate E ; and the resulting i ntegro-di fferential form 
of the Boltzmann transport equation is given as 



V • 5 <t> (r,E,^) + p<j>(r,E,^) = S(r,E,^) 



Z s (r,E'+E, ft' -ft) <J>(r,E', $') dE 1 d 



E' fi 



( 2 . 10 ) 



1 0 



In the above equation, the quantity p is the linear attenu- 
ation coefficient and the angular flux <J> ( r , E ,ti) is defined 
so that 



$(r,E,6) dE d^ (2.11 ) 

gives the number of photons at r, with energy E in the 
range dE, going in the direction t within the element of 
solid angle d:*, which cross in unit time a unit differential 
element of area whose normal is in the direction The 

various terms in equation (2.10) are defined as follows: 

V • 6 d> is the loss rate due to leakage from the phase 
space cell (dr, dE, d^), 

p < p is the loss rate due to collisions in the phase space 
cell which either absorb the photons or change their energy 
and direction, 

S is the independent generating source of photons, and 

r 

Z (?, E'-E, <J>(r, E'n) dE' dQ' (2.12) 

• E' h 

represents the rate at which photons are scattered into the 
phase space cell. The transport kernel, Z ( r , E 1 -*-E 1 ) , 

represents the differential macroscopic cross section for 
scattering from the direction 5' to H and from the energy 
E’ to E per unit solid angle and energy range. For Compton 
scattered photons, the transport kernel is given by the 
well-known Klein-Nishina formula. This formula takes a 



particularly simple form if the energy dependence is ex- 
pressed in terms of photon wavelength in Compton units, 

m r 2 

C 511 

A(Compton units) = = £ ( Me V ) ’ (2.13) 

In these units, the transport kernel is 

2 

2 s (r,A ’-X.ft'-S) = £ r e 2 {^-) [^r + + 2 ( A 1 - x) 

+ (A 1 -A) 2 ] 6(1 + A’ - A - 3 ) , (2.14) 

where n is the electron density and r g is the classical 
radius of an electron (2.817 x 10 cm). The Dirac delta 
function 6 accounts for the fact that wavelength (energy) 
and angle changes are related by the Compton shift equation 

A* - A = 1 - 3' • ft • (2.15) 

In shielding problems where absorbed energy density 
and biological dose are the quantities of interest, it is 
convenient to work with the angular energy flux which is 
defined as 



I (r , E ,fi) = E • <ji(r ,E,$) . (2.16) 



* One unit of Compton wavelength = h/m C = .02425°A, 
where h is Planck's constant, m Q is the rest mass of an 
electron, and C is the velocity °of light. 



This preference arises from the fact that absorbed energy 
density and biological dose are more nearly proportional 
to energy flux than particle flux. The transport equation 
is written in terms of I(f,A,^) by multiplying equation 
(2.10) by E and noting that 



and 



then, 



E • S(r,E,5) = S(r,A,$) 

I(r,E,$) = I (r.A.ft) 

, 2 

Z $ (r, E'+E,sW3) = A.ft'-H) ; 



(2.17) 

(2.18) 

(2.19) 

( 2 . 20 ) 



V • ^ I(r,A,fi) + u I(r,A,$1) = S(r,A,Q) 



A 



I (r, A'-A^'-ft) I (r , A ' ) d A ' d*'.(2.21) 



A direct, analytical solution of the Boltzmann trans- 
port equation is feasible only for a very small number of 
simple and highly idealistic problems. There are, however, 
many methods available which provide approximate solutions 
to the transport equation, several of which are described 
below. In general, methods utilizing numerical techniques 
with high speed computers provide more accurate solutions 
than methods which employ simplifying assumptions to alter 
the transport equation such that an analytical solution is 
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obtainable. An example of the latter type is the elementary 
one-speed diffusion equation of reactor physics which can 
be derived from the transport equation by invoking the 
following assumptions: 

1. Scattering is isotropic in the laboratory co- 
ordinate system. 

2. There are no strong absorbers in the medium 
i • e . , Y. a << z s . 

3. The region of interest is two or more mean free 
paths from any intense sources or boundaries. 

The existence of these limiting assumptions (especially the 
first assumption) is a clear indication that diffusion 
theory is inadequate for problems involving photon trans- 
port because the scattering of photons is not isotropic at 
any energy. 

2.1.1. The Method of Successive Scattering 

A relatively simple and straightforward approach for 
calculating the contribution of scattered photons to the 
penetration of gamma-rays is found in the method of suc- 
cessive scatterings. ^ ^ This method is a semianalytical 
technique which considers the total photon flux to be the 
sum of the unscattered flux, the once scattered flux, the 
twice scattered flux, etc. In notational form 

oo 

<t>(x)= l <J> j ( x ) s 
i = 0 



( 2 . 22 ) 



where ^(x) is the flux at x which results from the i th 
scattering. As previously mentioned calculations for the 
unscattered flux are relatively simple to perform; however, 
this is not so for the scattered flux. In the successive 
scattering method, the unscattered flux provides the 
collision density for first scattering which is used as a 
source in calculating the once scattered flux. The process 
is then repeated with an integral recursion formula which 
gives a relationship between <j>. + .|(x) and ^.(x). In other 
words, the higher-order scattered fluxes are expressed as 
integral operators of the next lower-order scattered flux. 
After the first scattering, the mathematics becomes quite 
involved and there is a tendency to limit this type of 
calculation to only one or two orders of sea ttering . ^ ^ 
However, this simplistic approach yields unreliable results 
for deep penetration problems where many scatterings con- 
tribute to the flux and the once and twice scattered con- 
tributions are only a small part of the total flux. 

The method of successive scattering was first used by 
Peebles and Plesset^^’^] ^ Q ca i cu i a te the transmission 
of monodi recti ona 1 gamma-rays through finite slabs of iron 
and lead. The error in these calculations is estimated to 
be as high as 20% for slab thicknesses up to 20 mean free 
paths. A variation of the met.hod of successive scattering 
has been developed by Peebles^^ and elaborated upon by 

r 2i l 

Yarmush. It consists of considering the photon 



transmission through a thick slab to be composed of a suc- 
cession of transmissions through thin slabs which are taken 
thin enough to make the contributions of photons scattered 
three or more times negligible. Gamma-ray transmission 
in the thick slab is then obtained by successive matrix 
multiplication of the thin slab results. 

★ 

2.1.2. The Spherical Harmonics Method 
The method of spherical harmonics, as applied to the 
transport equation, consists of expressing the various an- 
gular dependent terms in a complete set of elementary 
functions such as a series of polynomials. In three- 
dimensional geometry', spherical harmonics are the logical 
choice; however, these reduce to Legendre polynomials for 
one-dimensional geometry with azimuthal symmetry . When 
the polynomial expansions are substituted into the trans- 
port equation, it separates into an infinite set of coupled 
differential equations which no longer involve the direc- 
tional variables. Practical methods of solution require 
that the infinite set be limited to a finite number of 
coupled equations, and this is achieved by truncating the 
polynomial expansions to some finite number of terms. For 



*Although not widely used today, this method embodies 
concepts which are used in other techniques described in 
this Chapter. It has therefore been included for peda- 
gogical reasons rather than its potential for problem sol- 
ving. 
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example, a P n (or nth order) approximation represents a 
solution in which the polynomial expansions are limited to 
n+1 terms. 

In general, the degree of anistropic scattering dic- 
tates the order of expansion required to yield accurate 
results. The P-j approximation is equivalent to diffusion 
theory which, as previously indicated, provides adequate 
results only when the scattered flux is nearly isotropic. 

It has been shown that a P ? approximation provides a vast 

[ 22 ] 

improvement over diffusion theory 1 - J and that the most 
severe anisotropy resulting from Compton scattered photons 
with energies less than 10 MeV can be adequately resolved 

r 2 3 1 r 2 4 1 

with a P ^ approximation. J Lanning L J has used the 
spherical harmonics method to investigate gamma-ray heating 
in multi-slab geometry. He successfully calculated the 
spatial distribution of the energy flux spectra out to ten 
mean free paths with a P^ approximation. 

2.1.3. The Moments Method 

T251 

The moments-me thod of Spencer and Fano L J is a semi- 
analytical technique which solves the Boltzmann transport 
equation in infinite geometry. With this method, as in 
the spherical harmonics method, the angular flux is expand- 
ed in terms of Legendre polynomials. *' In addition, spatial 
moments of the Legendre coefficients, I (x,E), which serve 
to separate energy and spatial dependence, are introduced 
by the definitions 



I £ ( x , E ) x m dx , 



1 7 



b m^ " 



m+ 1 



m ! 



+ oo 



( 2 . 23 ) 



where u Q is the linear attenuation coefficient at the source 
energy E q . In terms of these moments, the transport equa- 
tion reduces to a set of interlinked integral equations of 

r p c i 

the Volterra type L J with energy as the only independent 
variable. The spatial moments are straightforwardly eval- 
uated by numerically integrating these integral equations. 

It should be emphasized that the calculations to this point 
are rigorous because no assumptions have been made. The 
major source of error in this method results from recon- 
structing the spatial dependence of the flux with only a 
finite number of known moments. Usually six to eight 
properly chosen moments are sufficient to reconstruct the 
spatial flux with reasonable accuracy.* 

The method of moments is generally considered to be 
the technique of choice for problems involving photon trans- 
port in infinite media. Goldstein and Wilkins^^ have used 
the moments method to investigate in detail the penetration 
of monoenergetic gamma-rays in infinite homogeneous media, 
and they have used their results to develop an extensive 
compilation of energy flux spectra and buildup factors for 



*Th'e interested reader is referred to Reference [1] for 
a detailed description of two methods used to reconstruct 
the spatial flux The polynimial expansion method and the 
method of undetermined parameters. 
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infinite plane and point isotropic sources. Advantages of 

r 271 

the moments method, as stated by Golstein, J include: 

1. A detailed energy spectrum of the flux can be 
obtained. 

2. The required computing time for penetrations up 
to 20 mean free paths is relatively short. 

3. Most photon source configurations of interest 

are amenable to the method. (Spencer and 
r 2 p 2 9 1 

lamkin L ’ J have successfully employed the 
method to calculate the slant penetrations of 
photons in water and concrete.) 

4. Angular distributions of the flux may be obtained. 
(Certaine^^ used the moments method to calculate 
angular distributions from plane monoenergeti c 
sources . ) 

The most serious limitation of the moments method is 

that it is applicable only to infinite homogeneous media. 

Another disadvantage, less restrictive but still noteworthy, 

involves the difficulty in reconstructing and determining 

★ 

the flux behavior near the source 0 ess than one mean free 
path). Because of this, the method of moments is not a 
good technique for evaluating reflection coefficients (al- 
bedos). 



* T 3 1 1 

Berger L J has developed a method to alleviate the 
difficulty for point isotropic sources. 



2.1.4. The Method of Invariant Imbedding 

The invariant imbedding technique was originally devel- 
oped and used to calculate the reflection of diffuse light 
from a stellar atmosphere by the astrophysist Ambrazum- 
ian.^^’33] Bellman, Kalaba, and Wing^^ have extended 
the method into a general approach, and shown that it can 
be applied to a much broader class of problems which in- 
cludes neutron and gamma-ray transport. As applied to radi- 
ation, the invariant embedding formulation is not another 
method for solving the transport equation but rather a new 
fundamental approach to the transport of radiation in ma- 
terial media. 

Unlike the Boltzmann equation, the dependent variables 
of the invariant embedding formulation are reflection and 
transmission functions. Each of these parameters is de- 
fined by and thus satisfies a non-linear integro-differential 
equation which cannot be solved analytically. Hence, all 
of the solutions are numerical in nature. The reflection 
and transmission equations represent problems of the 
" i n i t i a 1 - va 1 ue " type. In contrast, the Boltzmann equation, 
as applied to photon transport in slab geometries, forms 
"boundary- val ue" type problems. In general, i ni ti al -va 1 ue 
problems are more amenable to numerical solution with a 
computer. As a consequence, the invariant embedding formu- 
lation has a computational advantage over the Boltzmann 
approach . 
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Mizuta and Shimizu^^’^^ have applied the method of 
invariant imbedding to the problem of photon transport in 
slabs and demonstrated that it is competitive in accuracy 
to other established methods such as Monte Carlo and the 
method of moments. Compared to stochastic methods, in- 
variant imbedding has the advantage that it requires less 
computing time. Compared to the moments method, it has 
the advantage of being applicable to slab geometries. The 
disadvantages of the method are that it is difficult to 
apply to other than slab geometry and that it is inefficient 
for thin shields. (Thin shields require almost as much com- 
puting time as thick shields.) 

2.1.5. The Discrete Ordinates Method 

G. C. Wick^^ in 1 943 was the first to suggest the use 

of discrete ordinates (or discrete directions) in transport 

r 3 s "1 

theory. In the years that followed, S. Chandrasekhar 1 J 
further developed the method and applied it to the study of 
radiation transport in stellar atmospheres. Early applica- 
tions of the discrete ordinates method were primarily 

limited to the transport of isotropically scattered radia- 

[391 

tion in slab geometries. Carlson, with some success, 

applied discrete ordinates to curvilinear geometries by 
using an " approximation." In the version, the angular 
dependence of the flux is approximated by a series of con- 
nected straight line segments; hence, the letter S signifies 
segments and the subscript N indicates the number of segments 
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used to construct the angular dependence. It was subse- 
quently recognized by Carlson and Lathrop^^^ that the 

approximation was only a special case of a more general 
formulation and this led to the development of the discrete 
ordinates method. 

The discrete ordinates method, as presently employed, 
is a means of numerically solving the Boltzmann transport 
equation with finite difference equations. The finite dif- 
ference equations are formulated as flow balances for dif- 
ferential phase space cells (dV dE d?). The more recent 
versions of the method accurately treat anisotropic scat- 
tering by using a Legendre series expansion to describe the 
scattering cross section, ( r , E 1 -»-E , . The fundamental 

assumption in the discrete ordinates method is that the 
inscattering integral can be approximated by a quadrature 
scheme which evaluates the integral as the sum of a dis- 
crete distribution. In one-dimensional geometries a Gauss- 

r 4 1 1 

Legendre quadrature formula may be used to provide the 

discrete distribution, with discrete direction cosines cor- 
responding to the Gaussian zeroes. 

The discrete ordinates method has been shown to be 

a versatile and powerful method for accurately solving the 

transport equation in both one-dimensional geometries 

(DTF IV, ANISN^^) and two-dimensional geometries 

(DOT^^). It has been successfully applied to neutron and 

[45 46 4 7 1 

gamma-ray deep penetration calculations 1 - ’ ’ as well 
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as to general core criticality problems. The results 

of these calculations and many others clearly indicate that 

the discrete ordinates method, as applied to shielding prob 

★ 

lems, has the following advantages: 

1. Depending somewhat on the sophistication desired, 
the calculations are easy to prepare. 

2. The method is not stochastic and flux errors 
at deep penetration are systematic rather than 
statistical . 

3. Secondary gammas may be calculated by the same 
method, either as a second calculation or simul- 
taneously with the neutrons. The gamma yield 
distribution may be dependent on the neutron 
energy at capture. 

4. The neutron energy range from highest fission ener 
gies to thermal, including upscatter, may be cal- 
culated by the same method. 

5. The one-dimensional calculations are much faster 
(in computer time) than similar Monte Carlo cal- 
culations. In two dimensions the type of problem 
and the desired answers determine whether or 
Monte Carlo is better. 



*These a 
F. R. Mynatt 



jljj^tages 



were compiled from a listing given by 
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The major disadvantage of the method is that it is not 
applicable to three-dimensional geometries. 

The discrete ordinates method is one of the methods 
used in this work to calculate the penetration of photons 
in two-layered slabs; therefore, a derivation is given in 
Appendix A which establishes the equivalence of the one- 
dimensional discrete ordinates equations with the transport 
equation and which also shows how the Gaussian quadrature is 
employed to evaluate the inscattering integral. 

2.1.6. Stochastic Methods 

Many physical processes including the emission and 
transport of radiation are probabilistic in nature. As a 
result, many types of problems in mathematics and physics 
can be solved by random sampling or stochastic techniques. 

The Monte Carlo method, as applied to radiation transport, 
represents a system of probabilistic solutions in which 
probabilities relating to microscopic emission and interaction 
processes are used to predict macroscopic phenomena such as 
flux and dose-rate. Simply stated, the Monte Carlo approach 
consists of constructing a series of analogue experiments 
(histories) which mathematically simulates the emission and 
transport of particles, and then analyzing the outcomes of 
a large number of histories to obtain relevant quantities. 

Each history (or experiment) includes the generation of a 
particle, its random walk through the medium, and its death 
which results when the particle is absorbed, exceeds some 
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physical energy or time boundary, or when its importance 
becomes insignificant. 

The method of random sampling was first applied to the 

solution of radiation transport problems by Ulam and Von 
T 50 5 1 1 

Neumann 1 - ’ J in the mid-1940's. In the early stages of 
development, the Monte Carlo method was based on direct 
simulation of the transport problem and was therefore found 
to be prohibitively costly and time consuming in calculat- 
ing deep penetrations with reasonable accuracy. In the 
years that followed, several techniques were developed to 
reduce the variance associated with the statistical behavior 
of a particle history. Some of these variance reducing (or 
biasing) techniques, as described by Kahn, ’ J include: 

1. Importance Sampling: The selection of events 

(source parameters, path lengths, etc.) or omis- 
sion of events (absorption) is based on an altered 
probability distribution which allows the more 
important events to be sampled more frequently. 

The exponential transform is one of the more widely 
used means of imoortance sampling because it is 
easier to alter selection of path length. Particle 
weighting is another widely used importance samp- 
ling technique in which absorption is not allowed 
but accounted for by reducing the weight of a 
scattered particle. 
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2. Use of Expected Value: A portion of the random 

walk process is replaced by its average value. 

3. Russian Roulette and Splitting: Unimportant his- 

tories ( i . e . , those which are unlikely to contribute 
to the answer) are terminated with some probability 
while important histories are multiplied to give 
several semi -independent estimates. 

In general, Monte Carlo is not the technique of choice 
for one-dimensional steady-state transport problems because 
other methods such as discrete ordinates are usually much 
faster (in computer time) and just as accurate. In two- 
dimensions the type of problem (deep penetration, irregular 
geometry, ducts) determines whether Monte Carlo or discrete 
ordinates is better; but, there is no competitor to Monte 
Carlo for obtaining rigorous solutions to three-dimensional 
or time-dependent transport problems. 

The Lattice Model Concept of Stochastic Transport was 

developed and successfully applied to neutron transport by 

T54l 

H. E. Hungerford. J It differs from conventional Monte 
Carlo methods in that the random walks of particles are 
constructed by sampling their motion at predetermi ned points 
using a direct simulation of the transport process. In 
the "Lattice Model Concept" space can be envisioned as being 
made up of a large number of small cubes whose sides are 
equal to a distance called the unit lattice distance. The 
method of sampling particle motion within this framework is 
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best described by H. E. Hungerford L J in the following 
excerpt: 



If we now restrict the locations of the sampling 
points to be only at the corners of the cubes, and 
nowhere else, we have already chosen our allowed 
directions of emission, namely, these directions 
represented by lines which can be drawn from any 
corner of the unit cube to any other, neglecting 
the repeats. There are altogether 26 allowed 
directions of motion, outward from any given source 
point (considered at the origin)--6 along the 6 
axial directions, 12 in-plane diagonal directions 
(the planes being defined by any 2 of 3 axes), and 
8 skew diagonal directions (to opposite corners of 
the cube). Actually, to show all these directions 
emanating from one point, one has to make a larger 
cube from 8 unit cubes and place the source at the 
central point common to all cubes. These 8 cubes 
correspond to the 8 octants of a sphere. The larger 
cube [may be] referred to as the representati ve cube. 

[551 

A. Razani L J applied the lattice concept to gamma-ray 
transport by developing a stochastic gamma-ray transport 
code called PUGT I. Another version of the code, PUGT II, 
which uses particle weighting will be used in this work to 
calculate the penetration of photons in two-layered slabs, 
and thereby provide independent verification of the discrete 
ordinates results; therefore, a derivation which relates 
the transport of photons in the lattice model to the inte- 
gral form of the transport equation is given in Appendix B. 
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3. THEORY OF PHOTON TRANSMISSION AND 
BUILDUP IN A TWO-LAYER SLAB 

3.1. Analytical Formulation of the 
Two- Medium Buildup Factor 

Given two slabs of different materials and a plane, 
monoenergeti c gamma - r ad i a t i on source incident on one of them, 
the problem under consideration is to develop a relatively 
simple expression which will predict the buildup and there- 
fore the transmission of radiation through the slabs, with- 
out having to resort to expensive and time-consuming com- 
puter calculations. The choice of the radiation quantity 

★ 

selected as the prime variable is the angular energy flux 
density, I(b^,E,^) in MeV per cm^ per second per steradian 
and per unit energy, where 

I(b 1 ,E,fi) - E • <j>(b 1 ,E,ft) 

and <j>(b-|,E,^) is the angular photon flux density. (The num- 
ber of photons of energy E at b-| passing through a unit of 
area in the direction Q , per unit time, per steradian, and 
per unit energy.) The analytical formulation to be 

*Quantities used in this development are defined in the 
Nomenclature preceding the text. 
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developed is based on the buildup factor concepts, as given 
in Chapter 2, and the geometrical configuration shown in 
Figure 3.1. More rigorously stated and referring to Figure 
3.1, the objective of the formulation to be developed is 
to determine the buildup factor which when used with Lam- 
bert's Law correctly predicts the energy flux density at 
the detector (point C) resulting from a plane, monoener- 
getic, monodi recti ona 1 photon source at point A. The two- 
medium buildup factor , as derived, will be a function of 
quantities which characterize the angular energy flux at 
the interface between the two slabs. As nearly as possible 
the formulation for two-layer buildups will follow that for 
one layer. The single-medium buildup factors previously 
defined contained two terms, one describing the scattered 
energy flux and, the other describing the unscattered ener- 
gy flux. The same scheme will be followed here. In addi- 
tion, a distinction will be made between photons which 
scatter in the first layer and those that scatter in the 
second layer. In order that the reader fully understand 
the terminology and import of this development, a series 
of terms used herein is defined. The superscript s is used 
to refer to the scattered radiation and the superscript u 
refers to the unscattered radiation. These superscripts 
may be used alone or in combination. Thus: 

1. I s is the scattered energy flux, 
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c 

DETECTOR 



Figure 3.1. Geometrical Configuration for a Two-Layer Slab. 
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2. I u is the unscattered energy flux, 

3. I su refers to the energy flux which has scattered 
in the first layer but not in the second layer, 

4. I us refers to the energy flux which doesn't scat- 
ter in the first layer, but scatters in the second 
1 ayer , and 

5. I ss refers to the energy flux which scatters in 
both layers. 

These and other terms will be more fully described as the 
detailed development unfolds. 

The scattered energy flux at the detector (point C) 
may now be written as the sum of three components 

^(b-j+b^) = I^Cb^+b^) + I us (b-|+b2) + I^tb-i+b^) — > 

cm sec 

(3.1) 

where I s u ( b -| + b 2 ) is the scattered energy flux at the de- 
tector arising from the scattered energy flux at the inter- 
face which doesn't scatter again in traversing the second 
layer, I us (b-|+b 2 ) is the scattered energy flux at the de- 
tector arising from the unscattered flux at the interface 
which scatters at least once in traversing the second layer, 
and I (b-|+b 2 ) is the scattered energy flux at the detector 
arising from the. scattered energy flux at the interface 
which scatters at least one more time in traversing the 
second layer. Here, as in the rest of this work, the 



i, 



31 



quantity b-j denotes the number of mean free paths in the 
first slab at the source energy E q and b^ is analogously 
defined for the second slab. That is, b, = p^(E Q )x^, and 
b 2 = R 2 ^ o ^ X 2 w ^ ere anc * ^2 are ^ e ^ lnear attenuation 
coefficients of the first and second media. 

The relationship between the energy flux I^^+b^) 
and the scattered angular energy flux at the interface, 

I S ( b ] , E , 5) , is given by 



I SU (b ] +b 2 ) = 



dE 



- Uo ( E ) p 

I (b^E.S) e c d — 2 



MeV 



cm -sec 



4 TT 



(3.2) 



In slab geometry the angular energy flux can be expressed 
as 



I S ( b 1 , E ,!^) = 



I (b 1 ,E,w) 

2 TT 



7 



MeV 



cm -sec-MeV-STER . 



(3.3) 



and the differential solid angle is given by 



d^ 



2 TT r d r 

w ^ 

P 



steradians , 



(3.4) 



where w = cose, and p and r are shown in Figure 3.1. Sub- 
stituting these relationships into equation (3.2) yields 



I SU (b 1 +b 2 ) 




J 0 



CO 



dE 



' o 



t s / h r u \ ~ u 2^ E ^ p w r dr 



MeV 

cm -sec 



(3.5) 
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The integration variable may be changed to w by using the 

2 3 

relations p = x^/w, r = tanO , and rdr = - x^ dw/w . Re- 
ferring again to Figure 3.1, one sees that w -*■ 0 as r -+ °°, 
and w -* 1 as r -*■ 0. Substitution of these variables and 

limits into equation (3.5) yields 
E 1 



I su (b 1 +b 2 ) 



dE 



I (b 1 ,E,w) e 



-b 2 (E)/w 



dw 



MeV 



cm -sec 



(3.6) 



where b 2 (E) = u 2 (E)x 2 . 

In proceeding further it will be necessary to expand 
the scattered, angular energy flux, I s ( b -j , E , w ) 4 in a series 
of Legendre polynomials. For a plane source, the Legendre 
series expansion takes the form 



I S ( b 1 , E , w ) = l I S £ (b 1 ,b 2 ,E) P £ ( W ) • (3.7) 

£=o 



It can be shown that beyond a certain number of terms i = m, 

there is a negligible contribution to the scattered angular 

energy flux; therefore, the series is truncated at m where 

★ 

m indicates the order of the expansion. The terms P f (w) 
are the Legendre polynomials of order £, and the expansion 
coefficients are given by 



lS £ (b r b 2,E) 



1 

I S ( b, , E , w ) P (w) dw — A-- eV • (3.8) 

, 1 crn -sec-MeV 



*It is anticipated that a fifth order expansion will be 
sufficient to describe I s (b-|,E,w) to very good accuracy. 
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The physical significance of the expansion coefficients can 
be ascertained by noting that for z = 0 and z = 1 , these quan- 
tities represent the scattered energy flux spectrum and the 
scattered energy current spectrum at the interface between the 
two layers. The higher % terms define moments of these. In 
addition, the scattered energy flux at the interface is relat- 
ed to the zeroeth order Legendre expansion coefficient by 

fE 



I S ( b • 



b 2 ) - 



I^(b 1 ,b 2 ,E) dE 



MeV 



cm -sec 



(3.9) 



Substituting the Legendre expansion for I s ( b ^ , E , w ) into 
equation (3.6) gives 



m 



fE 



I SU ( b^ +b 2 ) - l 



1-0 



dE 



I ^ ( b i ,b 2 ,E) P £ (w) e 
MeV 



-b 2 (E)/w 



dw 



T 



. (3.10) 



cm -sec 



For convenience in later analyzing and tabulating, it is de- 
sirable to consider a two-medium parameter which is a slowly 
varying function of the thickness of the first slab. This 
quantity results when both sides of equation (3.10) are di- 
vided by the scattered energy flux at the interface as de- 



V b l ,b 2 ,E) , . -b 2 (E)/w 

— — — p £ ( w ) e dw 



fined in equation (3.9). 


Thus, 


m 


fE f 


T S U / k _l k \ v 2&+1 

>„ (», + b 2 ) - l 2 


0 

dE 


Z= o 


0 



o 1 ( b l ,b 2^ 



(3.11 ) 
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where 1 u ( b ^ + b£ ) is nothing more than I s u ( b -j + b^ ) normalized 

2 

to a scattered energy flux at the interface of IMeV/cm -sec. 

For convenience in carrying out the integration with respect 

to w, let b = b 2 ( E ) . Note that b / b£(E ). Also, let us 

define a set of quantities 
[ 1 



M fb) 
l ' 



P (w) e b/w dw for z = 1,2,. ..,m. (3.12) 

0 



Substitution of M ( b ) into equation (3.11) yields 



m 



C'VM - I M 11 



2 = o 



E o I ^ ( b 1 ,b 2 ,E) 
o I$ ( b l’ b 2 ) 



M„(b) dE. 



(3.13) 



The evaluation of M^(b) is accomplished by comparing it to 



the exponential integrals which are defined as 

I 

- V 



E„(b) = b 



n-1 



; 



d v 



n 



(3.14) 



o 



The relationship between M 0 (b) and the exponential integrals 
is made clearer by substituting v = b/w into equation (3.12) 
so that 



V b) = 



P fbv !_" V b dv 
i v v ‘ 



T 



(3.15) 



Specifically for i = o, P ( b / v ) = 1, and 



M 0 ( b ) . - E 2 (b) , 



(3.16) 
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for £ = 1, P^(b/v) = b/v, and 

M ] (b) = E 3 (b), (3.17) 

for £ = 2, P 2 (b/v) = ^-[3(b/v) 2 - 1], and 

M 2 (b) = 1[3 E 4 (b) - E 2 ( b ) J , (3.18) 

etc . 

Since P ( w ) is a polynomial in w , M(b) can be expressed 

£ £ 

in shorthand form as 

M £ (b) = P f [E j+2 (b)] • (3.19) 

Here the somewhat strange notation signifies that whenever 
w J occurs in P^(w), it is to be replaced by E. +2 (b). For 
example, when £ = 3, 



P 3 (w) 


= ^[5 w 3 - 3 w] , 


(3.20) 


M 3 (b) 


= 1[5 E 5 (b) - 3 E 3 ( b ) ] . 


(3.21) 



With M^(b) now expressed as an exponential integral 
function of energy, the final step in evaluating I^ u (b-|+b 2 ), 
as indicated in equation (3.13), involves integration with 
respect to energy. In performing this integration it is 
necessary to specify the expansion coefficients, I ^ ( b-j , b 2 , E ) ; 
therefore, the energy flux I^ u (b-|+b 2 ) will be calculated in 
the next chapter after the expansion coefficients are de- 



termi ned . 
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The scattered energy flux at the detector which was 
previously expressed as 



I i (b 1 +b 2 ) = I su (b ] +b 2 ) + I Ui (b ] +b 2 ) + I"(b 1 +b 2 ) — J 



us 



s s 



MeV 



cm -sec 



( 3 . 22 ) 



can be displayed in a more convenient form by defining a 
scattered energy buildup factor, B^.(b-|,b 2 ), as the variable 
which gives the buildup of the scattered energy flux at the 
interface as it traverses the second medium. That is. 



B| ( b 1 ,b 2 ) = 



I su (b 1 +b 2 ) + I ss (b ] +b 2 ) 
I SU (b 1+ b 2 ) 



1 + 



I SS (b ] +b 2 ) 

I SU (b ] +b 2 ) 



( 3 . 23 ) 



With the scattered energy buildup factor defined in this 
manner, it follows that 

s / ,_ . U x _ n s / l. l x T s / 1 . \ ,su /L x . tUS/. .. x MeV 



I i (b 1 +b 2 ) = B|(b 1 ,b 2 ) r(b ] ,b 2 ) I* u (b ] +b 2 ) + I Ub (b 1 +b 2 ) — f 



cm -sec 



( 3 . 24 ) 



us 



Further simplification is achieved by expressing I in terms 
of the previously defined single-medium energy buildup fac- 
tor, that is 

( b i +b 9 ) 



I us (b,+b 2 ) = E 0 e ' 1 2 [B E (b 2 ) - 1] — 



( 3 . 25 ) 



cm -sec 
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where 



B E (b 2 ) = 1 



+ 



I us ( b 1 +b 2 ) 
I U (b 1 +b 2 ) 






and 

I U (b 1 +b 2 ) 



- (b, +b ? ) 

E o e 



MeV 

~ 7 

cm -sec 



(3.26) 



(3.27) 



Substituting this expression for I us into equation (3.24) 
yields 



I s (b ] +b 2 ) = B E ( b -j , b 2 ) I s ( b n , b 2 ) I^ u (b 1 + b 2 ) 



+ E e 
0 



-(b,+b 2 ) 



[B E (b, ) - 1] — £ eV 



(3.28) 



cm -sec 



In keeping with the definitions for s i n g 1 e -med i urn 
buildup factors, the two-medium energy buildup factor, 
E E^l + * 3 2^’ 1S defined as the ratio of total energy flux to 
unscattered energy flux. Specifically, in formula notation 



B E (b 1 +b 2 ) 



I U (b ] +b 2 ) + I S (b ] +b 2 ) 
I U ( b 1 +b 2 ) 



I S (b ] +b 2 ) 

-(b, + bj 
E e 1 c 



0 

Substituting the expression for I s (b-|+b 2 ), as 
equation (3.28), into equation (3.29) yields 



given 



(3.29) 



i n 
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B E (b ] +b 2 ) = B E (b 2 ) 



Br(b,,b ? ) I s (b, ,b ) I^ u (b +b ) 

+ — >■ . ■ . - ! — (3.30) 



E o e 



-(b ] +b 2 ) 



Equation (3.30) can be further simplified by noting that the 
energy buildup factor at the interface between the two layers 
is given by 



B e ( b.| , b 2 ) = 1 + 



I S (b ] ,b 2 ) 



-b 



E o e 



1 



(3.31 ) 



where I s (b-j,b 2 ) is the scattered energy flux at the inter- 
face. Thus, 



B^(b ,b )[B (b, ,b )-l] I p U ( b-i +b ? ) 
B E (b 1+t>2 ) - B E (b 2 ) ♦ -i — 1 - 2 Lp(-b 2 ) ° 3 ^.(3.32) 



This is the formulation which will be used to evaluate energy 
buildup factors for two-material combinations, other than 
those of water, iron, and lead which are calculated in the 
next chapter. It will be shown in the discussion of results 
that the two-medium parameters B E (b^,b 2 ) and I^ u (b^,b 2 ) are 
slowly varying functions of b^ and smoothly varying functions 
of b 2 which can be accurately evaluated for many two-material 
combinations by simple linear interpolation over atomic num- 
ber . 

The energy buildup factor at the interface between the 
two layers, B E (b^,b 2 ), needs further clarification so that 
it can be evaluated from available single- medium parameters 
such as buildup factors and albedos. The scattered energy 
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flux at the interface is composed of the following fluxes: 

I ^ ( b i ) - The scattered energy flux which results 

at the interface when the second medium 
is replaced by a vacuum. 

- The scattered energy flux at the inter- 
face which is reflected from the second 
1 ayer . 

In notational form, 

I S (b 19 b 2 ) = I ^ ( b -j ) + lj(b 1 ,b 2 ) — ip • (3.33) 

cm- sec 

Now, the scattered energy flux which is reflected can be 
expressed in terms of flux albedos as 

lj(b 1 ,b 2 ) = I U (b-j ) a|J(b 2 ) + I^b^ a^(b 1 ,b 2 ) — 2 ~ — >(3.34) 

where the unscattered energy flux albedo is defined as 
the ratio of reflected energy flux to incident, unscattered 
energy flux; and the scattered energy flux albedo is de- 
fined as the ratio of reflected energy flux to incident, 
scattered energy flux. Since the unscattered energy flux 
is both monoenergeti c and monod i recti ona 1 , is equivalent 

to the energy flux albedos for normally incident sources 

r 5 1 

found in the open literature. J However, the scattered 
energy flux albedo is not dealt with in the open literature; 
therefore, it will be calculated and tabulated in this 



work . 
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Substituting equation (3.34) into equation (3.33) and 
dividing by the unscattered energy flux I u ( b ^ ) yields 

B E (b ] ,b 2 ) = B E ( b i ) + a^(b 2 ) + [B E ( ) - 1 ] a* ( b ] , b 2 ) , (3.35) 

where B^(b^) is the single-medium energy buildup factor for a 
normally incident source in slab geometry which has been 
evaluated in Reference [5]. It will be shown in the dis- 
cussion of results that the interface energy buildup fac- 
tor can be accurately evaluated for a large number of two- 
material combinations by using equation (3.35). 
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4. CALCULATIONS 
4.1. Genera 1 

The transport of monoenergeti c gamma-rays normally 
incident on two-layered slabs of water, iron, and lead has 
been investigated by using both discrete ordinates and sto- 
chastic methods to solve the Boltzmann transport equation 
in slab geometry. Source energies of 1, 3, and 6 MeV were 
chosen so that the two-medium buildup factor formulation, 
as derived in the previous chapter, could later be verified 
for the primary photon interaction processes: The photoelec- 

tric effect, Compton scattering, and pair production. The 
production of annihilation radiation has also been included 
in these calculations, but other secondary interaction pro- 
cesses such as fluorescence radiation and bremsstrahl ung 
have not. The consequences of not dealing with these photon 
interactions in the calculations will be considered and 
accounted for in the discussion of results. 

The discrete ordinates calculations were carried out 

with the Westinghouse version of a transport code called 

ANISN,^^] and the stochastic calculations were carried out 

r 5 5 1 

with PUGT II, J a gamma-ray transport code developed at 
Purdue University. On the average, the PUGT II stochastic 
calculations required much more computer running time than 
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the ANISN discrete ordinates calculations, especially for 

★ 

slab configurations of more than three mean free paths. 
Therefore, the stochastic calculations were limited to two- 
layered slab configurations of up to five mean free paths 
and primarily used to verify the accuracy of the discrete 
ordinates results. One of the major goals of this work 
will be to establish definitive accuracy limits for the two- 
layered buildup factor results. 

4.2. Stochastic Calculations 
Over 150 independent stochastic calculations have been 
carried out with the transport code PUGT II to evaluate the 
transmission and buildup of monoenergeti c photons normally 
incident on two-layered slabs of water, iron, and lead. 

PUGT II (Purdue University Gamma-Ray Transport) is a sto- 
chastic computer code which is based on the "Lattice Model 
Concept" as discussed in Chapter 2 and Appendix B. In brief, 
the code generates source gamma-rays and runs analogue ex- 
periments in which the photons are followed through a ma- 
terial medium until they exceed some physical, energy, or 
time boundary. The record of each photon from birth to 
death is called its history, and quantities of interest such 
as fluxes, currents, and albedos are obtained by summing 
the results of many thousands of histories. Reference (55) 
should be consulted for a more detailed description of PUGT n. 

*A typical ANISN run for slab configurations up to 13 
mean free paths required approximately 250 seconds of time on 
the CDC 6500 computer at Purdue University while a PUGT II 
run for a slab configuration of 4 mean free paths required 
1500 seconds of computer time. 
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The results of these stochastic calculations for two- 
layered slabs of water, iron, and lead will be primarily 
used to determine the accuracy of the discrete ordinates 
calculations; therefore, it is imperative to estimate and 
limit the inherent statistical inaccuracies. For statistical 
purposes, the transmission of photons through a slab can be 
considered to be a series of independent Bernoulli trials 
where M, the fractional number of photons transmitted through 
the slab in N histories (trials), represents the probability 
that a photon will be transmitted through the slab on any 
given history. The standard deviation for this statistical 
quantity is then estimated by 



C _ r M(l-M) -, 1/2 

^ - L N J 



(4.1) 



and the fractional deviation which indicates the percentage 
of deviation from the mean value of the fractional number 
of photons transmitted is given as 

F = 100 ~ % (4.2) 

It is believed that a fractional deviation limit of 2% will 
provide a sound basis for determining the accuracy of the 
discrete ordinates calculations; therefore, the number of 
histories for each of the stochastic problems in this work 
was chosen large enough to limit the fractional deviation to 
within 2%. On the average 10,000, 20,000, and 30,000 

*The fractional deviation was less than 1.5% for problems 
in which the second material was not lead. 



44 



histories were required for slab configurations totaling 2, 

3, and 4 mean free paths respectively. 

It should be noted that the statistical model, as 
presented above, does not account for the effects of biasing. 
Particle weighting is the only biasing technique used in the 
PUGT II code. This weighting procedure reduces the sta- 
tistical variance by analytically accounting for photoelectric 
absorption. However, the resulting decrease in the frac- 
tional deviation cannot be quantitatively determined. 

4.3. Discrete Ordinates Calculations 

Over 900 independent discrete ordinates calculations 
were carried out to calculate, as accurately as possible, 
the buildup and transmission of monoenergeti c gamma-rays 
normally incident on two- layered slabs. The bulk of these 
calculations dealt with two-layered slabs of water, iron 
and lead, source energies of 1, 3, and 6 MeV, and shield 
thicknesses of up to 5 mean free paths of the first layer 
followed by up to 8 mean free paths of the second layer. 

The materials water, iron, and lead- were chosen because they 
are commonly used together in pair combinations, and they 
represent the light, intermediate, and heavy elements re- 
spectively. The source energies of 1 , 3, and 6 MeV were 
chosen in order to account for buildup and transmission char- 
acteristics resulting from all of the primary gamma-ray 
interaction processes: The photoelectric effect, pair pro- 

duction, and Compton scattering. 
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The discrete ordinates calculations in this work were 
carried out with the one-dimensional, discrete-ordinates 
code ANISN. This computer program is written in 

FORTRAN IV and solves the one-dimensional Boltzmann transport 
equation in slab, cylindrical or spherical geometries. 

Theoretical considerations, as discussed in Appendix A, 
indicate that solution accuracy is highly dependent on the 
following parameters: 

1. Order of scattering (P r ). 

2. Order of angular quadrature (S^). 

3. Mesh interval spacing ( r). 

4. Cross section set used. 

The order of scattering (P„) represents the number of terms 
in the truncated Legendre series used to approximate the 
inscattering function. In general, the degree of anisotropic 
scattering dictates the order of scattering required to give 
accurate results. A (Pg) expansion was selected for the 
ANISN calculations in this work. 

The quadrature data sets used in the ANISN code con- 
sist of direction cosines and associated weights which spec- 
ify the discrete directions over which the inscattering 
function is evaluated. Symmetric sets corresponding to the 

Gaussian zeroes are given for various orders of angular 

T 5 7 1 

q ua d ra t ure . i n the ANISN User's Manual. J In general, for 
anisotropic scattering, the order of angular quadrature 
(S^) should be at least twice as large as the order of 



